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Abstract 

To understand the role of isospin effects on fragmentation due to the collisions of charge 
asymmetric nuclei, we have performed a complete systematical study using isospin dependent 
quantum molecular dynamics model. Here simulations have been carried out for 124 X n + 124 X n , 
where n varies from 47 to 59 and for Y m + 40 Y m , where m varies from 14 to 23. Our study 
shows that isospin dependent cross-section shows its influence on fragmentation in the collision of 
neutron rich nuclei. 
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I. INTRODUCTION 

Heavy-ion collisions have been extensively studied over the last decades. The behavior of 
nuclear matter under the extreme conditions of temperature, density, angular momentum 
etc., is a very important aspect of heavy- ion physics. Multifragmentation is one of the 
extensively studied field at intermediate energies. One of the major ingredient in heavy ion 
collisions is the symmetry energy, whose form and strength is one of the hot topic these 
days jlj. This quantity vanishes at a certain incident energy. Finite nuclei studies predict 
values for the symmetry energy at saturation of the order of 30-35 MeV. In heavy ion 
collisions highly compressed matter can be formed for short time scales, thus the study of 
such a dynamical process can provide useful information on the high density dependence 
of symmetry energy. Even at low incident energies which belong to even smaller baryonic 
densities, the isospin dependence of the mean field potential was shown to yield same result 
obtained with potentials that has no isospin dependences. These results are in similar lines 
and it also indicates that even binary phenomena like fission will also be insensitive towards 
isospin dependence of the dynamics [2j. Recently theoretical studies on the high density 
symmetry energy have been started by investigating heavy ion collisions of asymmetric 
systems 3j, |4|. Comparisons of collisions of neutron-rich to that of neutron-deficient 



systems provide a means of probing the asymmetry term experimentally 



-ti|. The 



experimental analysis of the isospin effects on fragment production has yielded several 
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interesting observations: Dempsey et al. 8] in their investigation of 124 ' 136 Xe -f 112 > 124 Sn 
at 55 MeV/nucleon found that multiplicity of IMF's increases with the neutron excess 
of the system. A more comprehensive study was carried out by Buyukcizmeci et al. 9] 
showed that symmetry energy of the hot fragments produced in the statistical freeze-out 
is very important for isotope distributions, but its influence is not very large on the mean 
fragment mass distributions. Symmetry energy effect on isotope distributions can survive 
after secondary de-excitation. Moreover Schmidt et al. [10[ in their investigation on 
the analysis of LCP's production and isospin dependence of l2A Sn + 64 Ni, 124 Sn + 58 Ni, 
12A Sn + 27 Al at 35 MeV/nucleon and 25 MeV/nucleon collisions found that isospin effects 
were demonstrated in the observables, such as the angular distribution of light particles 
emitted in central collisions at 35 MeV/nucleon and LCP's emission. On the other hand 



Tsang et al. 11] in their investigation of 112 Sn+ 124 Sn, 124 S , n+ 112 Sn systems at an incident 
energy of E=50 MeV/nucleon showed the effects of isospin diffusion by investigating 
heavy-ion collisions with comparable diffusion and collision time scales. They showed that 
the isospin diffusion reflects driving forces arising from the asymmetry term of the EOS. 
With the passage of time, isospin degree of freedom in terms of symmetry energy and 
nucleon-nucleon cross section is found to affect the balance energy or energy of vanishing 



flow and related phenomenon in heavy-ion collisions 121 ] . J. liu et al., 13J studied the 
effect of Coulomb interaction and symmetry potential on the isospin fragmentation ratio 
(N/Z) gas /(N/Z)uq and nuclear stopping R. They showed that Coulomb interaction induces 
important isospin effects on both (N/Z) gas /(N/Z)n q and R. However, the isospin effects of 
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symmetry potential and Coulomb interaction on (N/Z) gas /(N/Z)u q and R are different. 
Our present study will shed light on isospin effects on multiplicity of fragments produced 
in the collision of charge asymmetric colliding nuclei. We present microscopic study of 
isospin dependent nucleon-nucleon cross section on charge asymmetric nuclear matter. In 
this paper our aim is two fold, one is to look for effect of density dependent symmetry 
energy on fragmentation and second is to look for influence of isospin dependent and isospin 
independent cross-section on fragmentation due to collision of charge asymmetric nuclei. 
This study is carried out within the framework of isospin-dependent quantum molecular 
dynamics model that is explained in section-II. The results are presented in section-Ill. We 
present summary in section-IV. 



II. ISOSPIN-DEPENDENT QUANTUM MOLECULAR DYNAMICS (IQMD) 
MODEL 

Theoretically many models have been developed to study the heavy ion collisions at in- 



termediate energies. One of them is quantum molecular dynamical model (QMD) 14], Il5| . 
which incorporates N-body correlations as well as nuclear EOS along with important quan- 
tum features like Pauli blocking and particle production. 

n past decade, several refinements and improvements were made over the original QMD 



14 , Il5| . The IQMD [16| model overcomes the difficulty as it not only describe the ground 
state properties of individual nuclei at initial time but also their time evolution. In order to 
explain experimental results in much better way and to describe the isospin effect appropri- 



ately, the original version of QMD model was improved which is known as isospin-dependent 
quantum molecular dynamics (IQMD) model. 



The isospin-dependent quantum molecular dynamics (IQMD) 16| model treats different 
charge states of nucleons, deltas and pions explicitly, as inherited from the VUU model. 
The IQMD model has been used successfully for the analysis of large number of observables 
from low to relativistic energies. The isospin degree of freedom enters into the calculations 
via both cross-sections and mean field. 

In this model,baryons are represented by Gaussian-shaped density distributions 

fi(f,P,t) = J- e-^-W) 2 ^ e-GHW^. (i) 

Nucleons are initialized in a sphere with radius R = 1.12 A 1 / 3 fm, in accordance with the 
liquid drop model. Each nucleon occupies a volume of h 3 , so that phase space is uniformly 
filled. The initial momenta are randomly chosen between and Fermi moment umfjop). The 
nucleons of target and projectile interact via two and three-body Skyrme forces and Yukawa 
potential. The isospin degree of freedom is treated explicitly by employing a symmetry 
potential and explicit Coulomb forces between protons of colliding target and projectile. 
This helps in achieving correct distribution of protons and neutrons within nucleus. 
The hadrons propagate using Hamilton equations of motion: 

dri d( H) dpi d(H) 



dt dpi dt dri 



(2) 
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with 



H)=(T)+(V 

-,2 



i 1 i j>i 

xfj{r,p,t)drdf*dpdtf. (3) 



The baryon-baryon potential V l \ in the above relation, reads as: 



ytj ( r _ Vg kyrme + Vy ukawa + V C J Qul + V^ m 

-f , — - 



expdf 1 -f\/n) ZiZje 2 
+ H- 



(\r' — r\/jj,) |f* — f] 

+t 6 -nTi5(n'-r J ). (4) 
f?o 

Here Zj and Zj denote the charges of i th and j th baryon, and T|, T| are their respective 
T 3 components (i.e. 1/2 for protons and -1/2 for neutrons). Meson potential consists of 

Coulomb interaction only. The parameters /i and t±, ,t$ are adjusted to the real part of 

the nucleonic optical potential. For the density dependence of nucleon optical potential, 
standard Skyrme-type parameterization is employed. The choice of equation of state (or 
compressibility) is still controversial one. Many studies advocate softer matter, whereas, 
much more believe the matter to be harder in nature. We shall use soft (S) equation of 
state that have compressibility of 200 MeV. 



The binary nucleon-nucleon collisions are included by employing the collision term of well 
known VUU-BUU equation. The binary collisions are done stochastically, in a similar way 
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as are done in all transport models. During the propagation, two nucleons are supposed to 
suffer a binary collision if the distance between their centroids 

h - Tj\ < y^-, °tot = a(y/s,type), (5) 

"type" denotes the ingoing collision partners (N-N, N-A, N-7r,..). In addition, Pauli blocking 
(of the final state) of baryons is taken into account by checking the phase space densities 
in the final states. The final phase space fractions Pi and Pi which are already occupied 
by other nucleons are determined for each of the scattering baryons. The collision is then 
blocked with probability 

Pbiock = 1-(1-Pi)(l-P 2 ). (6) 

The delta decays are checked in an analogous fashion with respect to the phase space of 
the resulting nucleons. 



III. RESULTS AND DISCUSSION 



To check the influence of density dependent symmetry energy on fragmentation we have 
simulated 50 Sn 124: + 50 Sn 124 and 50 Sn 107 + 50 Sn 12i reactions by using isospin dependent 
quantum molecular dynamics (IQMD) model at incident energy 600 MeV/nucleon for 



complete colliding geometry. The phase space generated using IQMD model has been 



analyzed using minimum spanning tree [MST] algorithm 16] and minimum spanning tree 



with momentum cut 171] [MSTP]. The results obtained are discussed as follows: 



In FigQ], shows multiplicity of free nucleons, LMF's and IMF's as a function of scaled 
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FIG. 1: Multiplicity of free nucleons, LMF's and IMF's as a function of scaled impact parameter, 
impact parameters for 50 Sn lM +50 Sn lM and 50 Sn 107 + 50 Sn 1M . Our findings are as follows: 

(1) As we move from central to peripheral collisions the number of free nucleons and LMF's 
decreases because the participation zone decreases which leads to the lower number of free 
nucleons and LMF's. But in case of IMF's, curve shows a "rise and fall" this is because for 
central collision the overlapping of participant and spectator zone is maximum so we get 
very small number of IMF's. For semi peripheral collisions the participant and spectator 
zone decreases so the production of IMF's increases and for peripheral collisions very small 
portion of target and projectile overlap so again few number of IMF's observed most of the 
fragments goes out as heavy mass fragments (HMF's). 

(2) Number of free nucleons, LMF's and IMF's in ^qSti 107 +50 SVz 124 produced are smaller 
as compare to 5 oSn 124 + 50 Sn 124 . This is because for the neutron rich system, heavy residue 
with low excitation energy will predominantly emit neutrons, a channel that is suppressed 
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in case of neutron poor nuclei. 
The equation 



E{ P ) = £( Po )(f) 7 (7) 
Po 

gives us the theoretical conjecture of how symmetry energy varies against density. 7, tells us 



the stiffness of the symmetry energy 18|. In Fig.l, the free nucleons, LMF's and IMF's for 
7 = and 7 = 0.66 both the curves clearly indicates the density dependence of symmetry 
energy. From fig.l, one can see that: 

(1) Small difference is observed in both curves in case of LMF's at a scaled impact parameter 
range from 0.0 to 0.4. This is because the density during the fragmentation is smaller than 
normal nuclear matter density. Hence the role of density dependent symmetry energy is 
negligible. 

(2) When we apply momentum cut in addition to space cut the number of free nucleons 
increases because at low impact parameter participant zone increases and large number of 
free nucleons produced on the other hand number of LMF's and IMF's decreases with MSTP 
cut. 

Although we have tried the simulation for two different parameterization of density depen- 
dent symmetry energy but influence of this is very small, because the density at which 
fragmentation take place is lower than normal nuclear matter density. Hence influence of 
symmetry energy on fragmentation is very small, which is in agreement with the observation 
ofref. [19 1 . 
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FIG. 2: Multiplicity of free nucleons and LMF's with N/Z. 



Now to check the role of different cross-sections on fragmentation for charge asymmetric 
colliding nuclei, we have chosen two set of reactions, one where mass of colliding nuclei is 
40, but charge varies from 14 to 23. For first set the chosen reactions are 40 X m + 40 X m , 
where 40 X m = ( 4( V 23 , 40 Sc 21 , 40 Ca 20 , 40 Ar 18 , 40 C7 17 , 40 5 16 , 40 P 15 and 40 St lA ) respectively. 
For second set we have chosen the reactions for which mass of colliding nuclei is 124, but 
charge varies from 47 to 59. Second set of reactions taken are l2i Y n + 124 Y n , where l24 Y n 



]2 *Ag 47 , 124 C(i 48 , 124 /r i49 , 124 ^ 50 , 124 / 53 , 124 Cs 55 , 124 Ba 5e and 124 Pr 59 ) respectively. All 



the simulations are carried out for b = 0.3 at 100 MeV/nucleon for symmetry energy corre- 
sponding to 7 = 0.66. Here we take three different nucleon-nucleon cross-section because at 
low energy cross-section have very large influence on fragment production. Moreover, they 



have small effect on fragment production for central collision, whereas 



is strongly influenced at semicentral (6=0.3 in this case) collisions 



ragment production 



20|. In FigfJl we have 
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displayed the multiplicity of free nucleons (FN's) and light mass fragments (LMF's) as a 
function of charge asymmetry (N/Z). One can clearly see the effect of different cross-section 
on the production of FN's and LMF's. It has been observed that: 

(1) If we fix a nn = a pp = a np = 55mb, then maximum production of FN's and LMF's takes 
place. Isospin effect can be clearly seen when we use isospin dependent cross-section cr iso 
(o"n P = 3cr ran = 3a pp ) and isospin independent cross-section a no i SO (a nn = <j pp = a np ). 
&noiso will reduce the cross-section and thus the number of collisions, hence lead to less pro- 
duction of FN's and LMF's but <j iso will enhance the number of collisions and hence the 
production of FN's and LMF's. Moreover, one can see that nearly constant difference in the 
production with Oi so and <J no i S0 . 

(2) Minimum production takes place for the case when N/Z = 1 i.e symmetric charge colli- 
sions. As we know that nuclei offer very interesting isospin situation where, the symmetry 
potential, Coulomb interaction and isospin dependent nucleon-nucleon collisions are simul- 
taneously present. The Coulomb interaction is an important asymmetry term which can 
bring an important isospin effect into the observable quantities in the intermediate energy 
heavy ion collision. 

(3) The symmetry energy term affects the LMF's more than that of free nucleons. The 



(N — Z) 2 plays a crucial role 21] . It has been studied that the isospin effects plays more 
important role in case of LMF's rather than free nucleons. 

From fig. 2, it is clear that the minimum of fragment production is achieved at N=Z in case 
of FN's because they are produced in collision dynamics. But for LMF's the fragment pro- 
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FIG. 3: Multiplicity of free nucleons, LCP's and IMF's with energy at fixed scaled impact parameter 
for 50 5n 124 +50 Sn 124 and 50 Sn 107 +50 Sn 124 . 

duction is nearly constant. 

Figj3l shows the variation of multiplicity of free nucleons, LCP's and IMF's with energy for 
central collision, for 50 Sn 124 + 50 Sn 124 and 50 Sn 107 + 50 Sn 124 for two different nucleon-nucleon 
cross-sections. It has been observed that multiplicity of free nucleons and LMF's increases 
with increase in energy. On the other hand, one can see a "rise and fall" in the multiplicity 
of IMF's; this behaviour is similar to the behaviour shown by Aladin group [22|. Moreover, 
number of free nucleons and LMF's produced is very large as compared to IMF's this is 
because for central collisions, interactions are violent so large number of free nucleons and 
LMF's produced. It is clear from the figure that slope of the curve is steeper in case of 
50 Sn 124 +50 Sn 124 than so^n 107 + 50 Sn 124 and this theoretical observation is in agreement 



with the experimental observation of Sfienti et al. 22]. This rise is due to the fact that in 



case of neutron rich system, heavy residues with low excitation energy will predominantly 
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FIG. 4: Multiplicity of IMF's as a function of Zbound- 
emit neutrons, a channel that is suppressed in case of neutron-poor nuclei. Here one can see 
the difference in the production of FN's, LMF's and IMF's due to different cross-sections. 
Since the proton number in both the cases is same but neutron number is different, so we 
expect some difference in the production. 

In FigJH we have shown IMF's as a function of Z bound . The quantity Z bound is defined 
as sum of all atomic charges Zi of all fragments with Z{ > 2. Here we observe that 
at semi peripheral collisions the multiplicity of IMF's shows a peak because most of 
the spectator source does not take part in collision and large number of IMF's are 
observed. In case of central collision the collisions are violent so there few number of 
IMF's observed and for peripheral collisions very small portion of target and projectile 
overlap so again few number of IMF's observed most of the fragments goes out in heavy 
mass fragments (HMF's). In this way we get a clear "rise and fall" in multifragmentation 
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emission. But the influence of <Ji SO and a noiso is negligible here because IMF's are produed 
from the participant zone. It is observed that IMF's shows the agreement with data at 
low impact parameters but fails at intermediate impact parameters due to no acess to filters. 

IV. SUMMARY 

By using isospin dependent quantum molecular dynamics model we have studied the 
role of isospin effects on fragmentation due to the collisions of charge asymmetric nuclei. 
Here calculations were carried out for 124 X n + 124 X n , where n varies from 47 to 59 and for 
40 Y m + 40 Y m , where m varies from 14 to 23. It has been observed that isospin dependent 
cross-section shows its influence on fragmentation in the collision of neutron rich nuclei and 
there is a constant difference in the production of FN's and LMF's with a iso and a noiso for 
charge asymmetric nuclei. 
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